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Abstract. We show that the attenuated geodesic ray transform on two 
dimensional simple surfaces is injective. Moreover we give a stability 
estimate and develop a reconstruction procedure. 



1. Introduction 

The geodesic ray transform, that is, the integration of a function along 
geodesies, arises as the linearization of the problem of determining a con- 
formal factor of a Riemannian metric on a compact Riemannian manifold 
with boundary from the boundary distance function. This is the boundary 
rigidity problem, see [28J for a recent review. The standard X-ray trans- 
form, where one integrates a function along straight lines, corresponds to 
the case of the Euclidean metric and is the basis of medical imaging tech- 
niques such as CT and PET. The case of integration along more general 
geodesies arises in geophysical imaging in determining the inner structure 
of the Earth since the speed of elastic waves generally increases with depth, 
thus curving the rays back to the Earth surface. It also arises in ultrasound 
imaging. Uniqueness and stability for the case of integration along geodesies 
on simple manifolds (see precise definition below) was shown by Mukhome- 
tov [19j in the two dimensional case. Explicit inversion formulas in the two 
dimensional case were given in [23] for the case of constant curvature, and 
in the general case Fredholm type inversion formulas were given. 

In this paper we consider the case of the attenuated geodesic ray transform 
in two dimensions that we proceed to define. 

Let (M,g) be a compact 2D Riemannian manifold with boundary. The 
geodesies going from dM into M can be parametrized by the set d+S(M) = 
{(x, £) G TM ; x G dM, |£| = 1, (£, u) < 0} where v is the outer unit normal 
vector to dM. For any (x, £) G d + S(M) we let t i-> j(t, x, £) be the geodesic 
starting from x in direction £. We assume that (M, g) is nontrapping, which 
means that the time t(x, £) when the geodesic j(t, x, £) exits M is finite for 
each (x,C) G d+S(M). 

If o G C°°(M) is the attenuation coefficient, consider the attenuated ray 
transform of a function / G C°°(M), 

l a f(x,0 = r (X ' /(7(^,£))expr f'a^^^ds] dL 



o 







Here (x,£) G d+S(M). 



i 
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A compact Riemannian manifold with boundary is said to be simple if 
given any two points in the boundary there is a unique minimizing geodesic 
joining the two points, and if the boundary is strictly convex. The notion of 
simplicity arises naturally in the context of the boundary rigidity problem 

m- 

Our first result shows that the attenuated ray transform on simple surfaces 
is injective for any attenuation coefficient. 

Theorem 1.1. Let (M, g) be a simple 2D manifold, and let a be any smooth 
complex function on M. If / is a smooth complex function on M such that 
I a f = 0, then / = 0. 

Moreover we will give stability estimates and a reconstruction procedure 
to recover / from its attenuated ray transform I a f. 

In the case where M = M 2 with the Euclidean metric, the corresponding 
injectivity result for the attenuated X-ray transform has been proved by 
different methods in Arbuzov, A. L. Bukhgeim and Kazantsev [2], Novikov 
[22], Natterer |20| . and Boman and Stromberg [5J. These methods also come 
with inversion formulas. If M is the unit disc in R 2 with Euclidean metric, 
a direct inversion formula was given by Kazantsev and A. A. Bukhgeim |12j . 
See Finch [9] and Kuchment [15] for surveys of these and other developments 
in Euclidean space. The Euclidean attenuated X-ray transform is the basis 
of the medical imaging modality SPECT. 

The attenuated geodesic ray transform arises in inverse transport prob- 
lems with attenuation [16] , [T7J , when the index of refraction is anisotropic 
and represented by a Riemannian metric. It also arises in geophysics where 
there is attenuation of the elastic waves. Rather unexpectedly, this trans- 
form also appeared in the recent works [TJ, [13] in the context of Calderon's 
inverse conductivity problem in anisotropic media. 

Although the attenuated ray transform is well understood in Euclidean 
space, much less is known about this transform on manifolds. Bal [4] proves 
injectivity and gives an inversion formula in the hyperbolic disc M 2 . Frigyik, 
Stefanov and Uhlmann [TU] prove injectivity when (M, g) is simple and g and 
a are real analytic, or close to real analytic. Sharafutdinov proves injectivity 
of the attenuated ray transform on manifolds with a condition involving a 
modified Jacobi equation in [25] , and the size and curvature of the manifold 
in [26]. Dos Santos Ferreira, Kenig, Salo, and Uhlmann [7] prove the analog 
of Theorem 11.11 on any simple manifold if ||ct||i°o(jw) is small. A similar 
result, with a slightly different smallness condition, also follows from the 
general stability theory of JJLOj - In these last results, the smallness condition 
arises since the methods involve a perturbation about the unattenuated case 
where a = 0. 

The inversion results for M 2 and H 2 , which assume no smallness condition 
for a, are based on complex analysis and holomorphic functions. We will give 
a geometric version of these complex analysis arguments for simple surfaces, 
thus establishing injectivity of the ray transform for arbitrary attenuation 
coefficients. One of the key tools will be the commutator formula for the 
geodesic vector field and angular Hilbert transform, established in |24] in 
the study of the boundary rigidity problem in two dimensions. 
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At this point, let us give some other results which follow from the methods 
presented here. The next theorem considers the attenuated ray transform 
for combinations of functions and 1-forms. If F(x,£) = f(x) + a(£) for some 
smooth function / and 1-form a, where a(£) = cij(x)£ 3 , the attenuated ray 
transform of F is defined by 



I a F(x,0= F(j{t,x,0,i(t,x,ti))exp / 0(7(5, x,0)ds 











dt 



where (x,£) € d+S(M). It is easy to see that this transform has nontrivial 
kernel since I a (ap + dp(£)) = for any p £ C°°(M) with p\qm = 0. The 
injectivity result, which also extends the corresponding result for functions, 
states that these are the only elements in the kernel. 

Theorem 1.2. Let (M,g) be a simple 2D manifold and let a € C°°(M) 
be a complex function. Suppose that / is a smooth function and a is a 
smooth 1-form on M, and let F{x, £) = fix) + aj(x)^ . If I a F = 0, then 
F = ap + dp(£) for some function p <G C°°(M) with p\qm = 0. 

Note in particular that if / = and a is nonvanishing, then any 1-form is 
uniquely determined by its attenuated ray transform. Results of this type 
were given in the unit disc in R 2 in |12| . for simple manifolds with ||a||x,°o 
small in [7], and for simple manifolds with g and a close to real analytic in 
|11| . For inversion formulas in M 2 see also [3], |21| . 

Once injectivity of I a is known, the general principle that the normal 
operator N a = (I a )*I a is an elliptic pseudodifferential operator and the 
arguments in [TO], [TJJ yield a stability result. To state this result properly, 
we use the solenoidal decomposition of a smooth 1-form a in M, 

a = a s + dp, 

where a s is solenoidal (meaning that 5a s = 0), and p\qm = 0. Here 5 
is the codifferential. This decomposition is uniquely determined by taking 
p = G{5a) where G is the inverse of the Dirichlet Laplacian on M. We also 
choose a simple manifold (M\,g) which is slightly larger than (M,g) and 
extend smooth functions and 1-forms in M by zero to M\. In this way iV a 
can be viewed as a pseudodifferential operator acting on functions in M\. 
The stability result is as follows. 

Theorem 1.3. Let (M,g) be a simple 2D manifold and let a G C°°(M) be 
a complex function. Suppose that / is a smooth function and a is a smooth 
1-form in M. Then 

||/ - aG(5a)\\ L 2 {M) + \\a s \\ L 2 {M) < C\\N a (f + a^)\\ 

Finally, we outline a reconstruction procedure to determine a function / 
from I a f. For simplicity we will assume that / is compactly supported in 
M mt and all quantities are real valued (the complex valued case is discussed 
in Section [6|). The reconstruction procedure consists of several steps, and 
we refer to Section [2] for a more precise explanation of the notations used in 
the result. 
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Theorem 1.4. Let (M,g) be a simple 2D manifold, and let a € C°°(M) be 
real valued. A real valued function / £ C£°(M mt ) can be determined from 
the knowledge of I a f using the following procedure: 

1. Define a function d on dS(M) by 

d(x = / 0. (x,Oea+s(M), 

1 0, otherwise. 

2. Find an odd holomorphic function w such that Jftfw = —a. 

3. Let (3 = (Id - iH)(e~ w d) on dS(M). 

4. Let v = /3 o V + u (/0) " 1(j41/3) in S*M, where A* = /3 - /? o ^ on 
d+S(M) and is the inverse of the geodesic ray transform in 
(M,g), in the sense that 

(I°)- 1 I ^ + a j e) = cl ) + a j e 

for a smooth function <p and a solenoidal 1-form a. 

5. Define rh = ^Re[(Id — iH)(e w v)] and u = rh — mo- 

6. Define q = (d — u)q o if) + ( u (^ Ju + au )- ) 0) an d let n = g + n. 

7. Let / = — (J^u + au)o. 

There are two nontrivial steps (Steps 2 and 4) in the above result: they 
require to find a holomorphic integrating factor w to the transport equa- 
tion (J£° + a)u = 0, and to invert the geodesic ray transform 1° with zero 
attenuation. Both these steps can be achieved in an explicit way if (M,g) 
has constant curvature, or if (M,g) is a small perturbation of a constant 
curvature manifold (see Section [6]). However, it is not clear how to carry 
out these steps explicitly in a general simple 2D manifold. 

It seems that even when M is a domain in M? with Euclidean metric, 
the reconstruction procedure does not reduce to a simple formula such as in 
|12j . |22| . It would be interesting to give a reconstruction procedure which 
would reduce to such a simple formula on constant curvature manifolds. 

The structure of the paper is as follows. Section [2] establishes notation 
and preliminaries related to geodesic flow, Hilbert transform, and functions 
which are holomorphic in the angular variable. In Section [3] we explain the 
strategy of the injectivity proof, starting with a simple inversion scheme 
based on holomorphic solutions of the transport equation, and discussing 
the modifications to this scheme required in the attenuated case. The first 
main step in the proof, the construction of holomorphic integrating factors, 
is achieved in Section [H using pseudo differential arguments. The second 
main step consists in proving that solutions of certain transport equations 
are necessarily holomorphic. This is done in Section [5l where also Theorems 
11.11 to 11.31 are proved. The final Section [6] gives a reconstruction procedure 
and proves Theorem 11.41 

Acknowledgements. M.S. is partly supported by the Academy of Fin- 
land. G.U. is supported in part by NSF and a Walker Family Endowed 
Professorship. 
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2. Preliminaries 

We refer to [23], [23], [25] for the following facts. Assume that (M,g) 
is a compact 2D Riemannian manifold with boundary dM. We will also 
assume that (M,g) is nontrapping and dM is strictly convex (see below). 
We denote the inner product on tangent vectors and other tensors by ( • , • } 
and the corresponding norm by | • |. 

2.1. Geodesies. We will mostly work on the unit sphere bundle given by 

SM = |J S x , S x = {(x,0 G TM; |£| = 1}. 

The manifold SM has boundary dS{M) = {(x,£) G SM; x G dM}. The 
outer unit normal vector of dM is denoted by u, and the sets of inner and 
outer vectors on dM are given by 

d±S(M) = {(x,£) eSM;xe dM, ±($,v) < 0}. 

If is a point in SM we denote by j(t,x,£) the geodesic on M 

satisfying 7(0, = x and j(Q,x,£) = £. The geodesic flow is the map 

ipt : SM^SM, <pt(x,$) = (j{t,x,£),i(t,x,Z)) 

if t is such that the right hand side is well defined. The nonnegative time 
when a geodesic j(-,x, £) exits M is denoted by t(x,£). The manifold 
(M,g) is said to be nontrapping if t(x, £) is finite for any (x,£) G SM. The 
boundary dM is said to be strictly convex if its second fundamental form is 
positive definite. 

Since (M, g) is nontrapping and has strictly convex boundary, the next 
result holds by [251 Section 4.1]. 

Lemma 2.1. r is continuous in SM and smooth in SM \ S(dM), and 
further the function r_ : dS(M) — > R defined by 

T ( T p\-{ %T(x,t), (x,0£d + S(M), 
{ ' U 1 -|r(x,-0, (*,£)€ d_S(M) 

is smooth. 

2.2. Scattering relation. The scattering relation a maps an inner unit 
vector (x, £) G d + S{M) to the outer vector £). Thus, a takes the 
starting point on the boundary and direction of a geodesic and gives out 
the endpoint and direction of that geodesic. It is possible to define a as a 
smooth map on all of dS(M) by 

Then a is a diffeomorphism dS(M) -> dS(M) and a 2 = Id. 

2.3. Geodesic vector field. The geodesic vector field J$? is the vector field 
on SM which acts on smooth functions u on SM by 

d 



THE ATTENUATED RAY TRANSFROM 



6 



We consider two boundary problems related to <2fc '. If F is a smooth function 
on SM, then the problem 

J4?u = —F in SM, u\d_s{M) = 

has the solution u = u F where 

u F (x,0 = [ ( ° F(ip t (x,0)dt. 



A + w(x,£) 



If to is a smooth function on d+S{M) then the problem 

fflu = in SM, u\q + s(m) = w 

has the solution u = given by 

Wip = w o a o tp (2-1) 

where tp is the end point map tp(x, £) = (p T r x £)(x, £), and a is the scattering 
relation. 

Since r is continuous on SM and smooth on SM \ S{dM), the same is 
true for u F and . It is a minor inconvenience that these functions are not 
smooth on SM in general. The space of those w for which is smooth in 
SM is denoted by 

C™(d+S(M)) = {w G C°°(d + S(M)) ; tity G C°°(5M)}. 

This space was characterized in [24J in terms of the operator A + of even 
continuation with respect to a, acting on w G C°°(d + S(M)) by 

w(x,0, {x,£)€d+S(M), 
w(a(x,0), (x,Oed„S(M). 

Lemma 2.2. C™(d + S(M)) = {w G C°°(d + ,S(M)) ; G C°°(95(M))}. 

As for u F , sometimes we can work under the extra assumption that F 
vanishes near dM in which case u F is smooth on SM. At other times, we 
can use the fact that the odd part of u F is smooth in SM provided that F 
is even. If u is a function on SM the even and odd parts are defined by 

u±(x,£) = -(u(x,£) ±u(x,-£)). 

Of course, u is called even (resp. odd) if u = u + (resp. u = 

Lemma 2.3. If F is an even smooth function on SM, then u F is a smooth 
function in SM and satisfies J4?u F = —F. 

Proof. The last statement follows since WJrffu = —J^R*u where R is the 
map R{x,£) = {x, —£)■ This implies that (J4?u) + = J4?U-. 

We will reduce the smoothness statement to Lemma 12.21 Let (M,g) be a 
nontrapping manifold with strictly convex boundary so that M C M int (this 
can be achieved by embedding (M, g) to a compact manifold {S, g) without 
boundary and by looking at a small neighborhood of M in 5). If r(x,£) is 
the exit time of geodesies in (M,g), we know that f is smooth in S(M mt ). 

Extend F as a smooth even function into SM, and define 

fi(x,0 = I ( ' € F(<p t (x,0)dt 
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where (ft is the geodesic flow in (M, g). Then u € C°°(SM) and Jtfu = —F 
in SM. 



u_ is 



Let uu = (u — u_)\q + s(m)- Since Jrff(u — u_) = in S(M ) and u 
continuous in SM, we obtain u — u^_ = w^. Thus, to show that is smooth 
in SM it is enough by Lemma I2.2I to prove that A + w is in C°°{dS{M)). A 
short computation, using that F is even, gives that for (x, £) G dS(M) 

A+w(x,0 = ~ F(<p t (x,Z))dt + - F(<p t (x,£))dt. 

1 JO 1 J2t-{x,C) 

We know that 2r_ is smooth in dS{M) by Lemma 12,14 hence also A + w is 
smooth. □ 

2.4. Hilbert transform. To discuss functions which are (anti)holomorphic 
in the angular variable, we introduce the fiberwise Hilbert transform which 
acts on smooth functions on SM by 

Hu(x, = ^J s 1 -J^^ U & dS ^), € SM. 

The integral is understood as a principal value. Here (£±)j = £jk£ k where e 
is the clockwise rotation by 90 degrees: 

1 
-1 



If Hq is the usual Hilbert transform on the unit circle S , and if F x is any 
orientation preserving isometry from S x onto S 1 (such a map is unique up 
to rotation on S 1 ), one has for fixed x 

H = F* x H (F x l )\ (2.2) 

The last identity allows to transfer standard properties of the Hilbert trans- 
form on the unit circle to the present setting (see also \17\ Section 8]). 

A crucial ingredient for our arguments is a commutator formula proved 
in [23], which gives a connection between the geodesic vector field and the 
fiberwise Hilbert transform. 

Proposition 2.4. If u is a smooth function on SM then 

[H, J^]u = J^ ± u + (J?lu)o- (2.3) 
Here uq is the average of u over the angular variable: 

M X ) = 7T I u(x,0dS x (£), xeM. 
27r Js x 

We have used the vector field Jtf± = (£j_) j Vj on SM. Here V is the 
horizontal derivative on SM [25]. In local coordinates it is given by 



We collect some further basic properties of the Hilbert transform 
These involve even and odd functions with respect to the angular variable. 
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Proposition 2.5. The Hilbert transform maps even (resp. odd) functions 
with respect to £ to even (resp. odd) functions. If u is a function on SM 
then Hu± = H±u where 

H + u(x, £) = -!- f J ,r, \ u(x,rj)dS x (ri), 
If 1 

H-u(x,£) = — j- ru(x,r])dS x {ri). 

2vr J Sx (£l,»7) 

Also, if « is a function on SM then {Hu)q = 0, and if it = u(x) then Hu = 0. 

2.5. Holomorphic functions. The arguments below will be based on the 
ability of finding (anti) holomorphic solutions to transport equations. Here, 
(anti)holomorphic refers to the angular variable. The precise definition uses 
the Hilbert transform and is as follows. 

Definition. A function u on SM is called holomorphic if 

(Id — iH)u = Uq. 
We say that u is antiholomorphic if 

(Id + iH)u = Uq. 

The next result, which follows by (|2.2p . will be used many times below. 

Lemma 2.6. The product of two (anti) holomorphic functions is (antiholo- 
morphic, and e w is (anti) holomorphic if w is (anti)holomorphic. 

As an example, and to obtain some intuition into the arguments below, 
we will discuss the above notions in the case where M is an open set in M? 
with Euclidean metric. Then SM = M x S 1 , and any function u on SM 
may be written as Fourier series 



u(x, e ld ) = u k{x)e l 



jke 

k=—oo 

Here Uk(x) are the Fourier coefficients 

n fc (x) = — e- ike u(x,e ie )d0. 
27T Jo 

Then the even and odd parts of u are obtained by just taking the even or 
odd Fourier coefficients, 

u + (x,e i6 )= Mx> ik6 , 

k even 

u.{x,e w ) = u k (x)e ikd . 

k odd 

Also, with the convention sgn(0) = 0, 

H(e ike ) = -sgn(k)ie ike . 
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Therefore 

oo 

(Id + iH)u = u (x) + 2 Mx)e ike , 

k=l 

-1 

(Id - iH)u = u (x) + 2 u k (x)e lke . 

k=~oo 

Now (ld±iH)u = uq means that the negative or positive Fourier coefficients 
vanish. Thus, u is holomorphic (antiholomorphic) if and only if for any x 
in M, u(x, ■ ) extends into a holomorphic (antiholomorphic) function in the 
unit disc. 



3. Strategy of proof 

Our proof of Theorem 11.11 reduces the attenuated ray transform to the 
analysis of solutions of a transport equation. Let be the geodesic vector 
field on the unit sphere bundle SM. Then I a f = u\q + sm where u satisfies 

(J^ + a)u = — / in SM, u\q_sm = 0. 

Holomorphic solutions of certain transport equations will be crucial in the 
proof. To explain why such solutions might be useful, we first discuss a 
simple scheme which would imply injectivity and which turns out to work 
in the unattenuated case. In the end of the section we outline the strategy 
for the attenuated case. 

First inversion scheme. Let a and / be real valued, and let u be the 
solution given above. Motivated by the earlier result [12] in M?, it turns out 
that injectivity of I a would be a consequence of the following idea: 

Produce a function u* , which is holomorphic (or antiholo- 
morphic) in the angular variable, such that (u*)o = and 

{J^ + a)u* = -f (3.1) 

and such that ii*|as(M) *s determined by u\g S ^ M y 

To see how the above statement could be used to invert the attenuated 
ray transform, it is enough to take the imaginary part of (|3,ip to obtain 

(je + a)(lmu*) = 0. 

bmce jj? + a = e u -J^e~ u - (recall that u a _ is smooth in SM by Lemma [2.3p 
this shows that e _n -Imu* is constant on geodesies, and therefore Imu* is 
determined by its boundary values on d+S{M). But u* is (anti)holomorphic 
with zero average so its real part is determined by the imaginary part, and 
we obtain that u* in SM is determined by u\q + s(m) = I a f ■ Then / can 
be reconstructed from I a f for instance by taking averages over the angular 
variable in (|3.ip : 

/ = -((jf + oK) . 

This would give an inversion formula for the attenuated ray transform. 
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The unattenuated case. The above scheme actually works in the case 
a = where no attenuation is present, and results in a similar inversion 
formula as in |23j. Let (M, g) be a simple surface, and assume for simplicity 
that / G C~(M int ) is real valued. We would like to recover / from the 
knowledge of the geodesic ray transform 1° f\d + s(M)- 
Let u be the solution of 

j^u = -f in SM, u| a _5(M) = 0. 

Then u = w and u\g + s(M) = I°f- To obtain a holomorphic solution with 
zero average, the first idea is to take 

u* = (Id + 1H>_. 

To see if u* solves the transport equation, we use the fact that Jifu- = — f 
and compute by (j2.3j) 

J^fu* = (Id + iH)J^u- - i[H, J>F}u- = -f - i(J?±u) . 

The last expression on the right was analyzed in [23\ Section 5]. 

Proposition 3.1. If (M,g) is simple then the operators 

wf = (jr ± ut) , 

W*f = (u^f)o 

have smooth integral kernels and extend as maps from L 2 (M) to C°°(M). 
Also, W* is the adjoint of W, and if (M,g) has constant curvature then 
W = W* = (this last statement is also true in the presence of conjugate 
points). 

It follows that 

JTu* = -/ - iWf. (3.2) 
Thus, if (M, g) has constant curvature, then Wf = and u* is the required 
holomorphic solution with zero average. Note that 

u *\dS(M) = ( Id + *#>-lc>S(Af) 

so u*|as(M) is indeed determined by u\g S ( M y The scheme above now gives 
an inversion formula for the geodesic ray transform 1°. 

In the case where (M,g) does not have constant curvature the quantity 
Wf may be nonzero, but one can still obtain a Fredholm type inversion 
formula as in [23 . The right hand side in (j3.2|) is complex so splitting into 
real and imaginary parts is not immediately useful. However, one can iterate 
once more and introduce the antiholomorphic odd function 

u** = (Id - iH)uL +iWf . (3.3) 

This satisfies by ([23]) and the fact that J^u f _ +iWf = -f - iW f 

Jtfu** = (Id - iH)Jfu f _ +iWf + i[H,J^]u f _ +iWf 

= -f- iWf + i(^±u f+iWf ) 

= -f- W 2 f. 

Now / + W 2 f is real and the inversion scheme above can be used to recover 
/ + W 2 f from I°f. Thus we have constructed / up to a Fredholm error. 
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The attenuated case. If (M,g) is a simple 2D manifold and a is a real 
attenuation coefficient, we will use a modification of this inversion scheme 
which still involves holomorphic solutions. We explain the idea for proving 
injectivity. Suppose that / G C^°(M mt ) is real valued and I a f = 0. Then 
the corresponding solution u of the transport equation satisfies 

(JT + a)u = -f in SM, u\ 9S m = 0. 

The first step in the proof is to find a holomorphic integrating factor: an odd 
holomorphic function w G C°°(SM) such that one has the operator identity 

e w J^e~ w = J^ + a. 

We derive a characterization for the existence of such w, and we will employ 
pseudodifferential arguments to show that one can always find a holomorphic 
integrating factor. Given such w, the function e~ w u satisfies 

J(e"\) = -e~ w f in SM, e- w u\ dSM = 0. 

The second main step is to show that any solution which satisfies such a 
transport equation with holomorphic right hand side and which vanishes on 
dSM is necessarily holomorphic. This uses the commutator formula f|2.3[) 
for H and Jff, and boils down to the injectivity of the unattenuated ray 
transform on 1-forms [T]. We obtain that e~ w u is a holomorphic function, 
and since e w is holomorphic then so is u. 

This argument shows that whenever I a f = 0, the solution u of the trans- 
port equation must be holomorphic. Since u is real it is also antiholomor- 
phic, which shows that u only depends on x. Thus u = uq, and the transport 
equation reads 

duo{0 + au = -f in SM, u \ dSM = 0. 
Evaluating at ±£ shows that duo = 0, and consequently no = and / = 0. 

4. Holomorphic integrating factors 

Let (M,g) be a simple surface and let a be a smooth function on M. 
We consider the problem of finding holomorphic integrating factors for the 
operator J%? + a. More precisely, we are looking for smooth holomorphic 
functions w on SM such that 

(Jf + a)v = e w J^{e~ w v) 

for all smooth functions v on SM. An equivalent condition is that w should 
satisfy 

Jtw = —a. 

The main result of this section shows that holomorphic integrating factors 
always exist. 

Proposition 4.1. Let (M,g) be a simple 2D manifold and let a be any 

smooth complex function on M. There exist a holomorphic w G C°°(SM) 
and an antiholomorphic w G C°°(SM), both odd functions, such that 
Jtffw = Jtffw = —a. 
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One smooth solution to the equation fflw = —a is given by w = u°_. We 
begin with the simple observation that on constant curvature manifolds, the 
projection of u a _ to holomorphic functions also satisfies this condition. 

Lemma 4.2. Consider the operators V, T acting on functions on M by 

Ta = (Id + iH)u°L , ta = (Id - iH)u a _ . 

Then T (resp. V) maps C°°(M) to odd functions in C co {SM) which are 
holomorphic (resp. antiholomorphic) in the angular variable. One has 

J^Ta = -a- iWa, J^fa = -a + iWa. 

If Wa = 0, then 

r J%?Ta = Jtf'Ta = —a. 

Proof. The function u°L is smooth by Lemma 12.31 hence Ta and Ta are odd 
and smooth. By the commutator formula (|2,3p 

J^Ta = (Id + iH)Jfu a _ - i[H, Jf}ut = -a - i(Jfiu a ) = -a - iWa. 

The computation for T is analogous. □ 

Corollary 4.3. If a G Ran(Id+iW / )nRan(Id— iW), then there exist smooth 
holomorphic w and antiholomorphic w such that Jftfw = fflw = —a. 

Proof. If a satisfies the given condition, there are b, b € C°°(M) such that 
a = b + iWb = b — iWb. Letting w = Tb and w = Tb, we have 

J? w = -b- iWb = -a 

and similarly for w. □ 

The corollary shows that on any manifold for which Id + iW and Id — iW 
are surjective, one can find holomorphic integrating factors. This includes 
constant curvature manifolds (since W = 0) and small perturbations of these 
(since W has small norm |14j). However, for general simple manifolds we 
do not know if Id ± iW are surjective and one needs to work harder. 

First we give a characterization of those a for which J^w = —a for some 
(anti)holomorphic function w. 

Lemma 4.4. Let a be a smooth complex function on M. The following 
conditions are equivalent: 

(1) There exists a holomorphic (resp. antiholomorphic) odd function w 
in C°°(SM) such that 

J4?w = —a. 

(2) There exists a function h € C^°(0+5'(M)) and an antiholomorphic 
(resp. holomorphic) even function b € C°°{SM) so that 

a = b + i(J%j_u b ) + i(J^fj_h^)o 

(resp. a = bo — i(J$?±u b )o — i(Jfih^)o). 
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Proof. Let first w be a holomorphic odd function with Jtffw = —a (the case 
of antiholomorphic functions is analogous). Then w = (Id + iH)w for some 
odd w E C°°(SM), in fact one may take w = 
Now by the commutator formula (|2.3p 

-a = J?w = Jf (Id + iH)w 

= (Id + iH)J^w - i[H, je]w 

= (Id + iH)M>w - i(Jfiw) . 

This shows that (Id+iH)Jifw only depends on x. Consequently the function 
b = —Jffw is antiholomorphic and even, and we have (Id + iH)J4?w = —bo. 
The equality J4?w = —b implies that J4?(w — u b _) = 0, hence w = u b ___ + 
for some h E C™(d+S(M)) by Lemma EES Thus 

a = b + i(Jfiu 6 )o + ip2j_/i,/,)o- 

We have proved (2). 

Conversely, assume a is of the form given in (2) with b antiholomorphic 
and even. Define w = u b _ + (h^)-. Then w is odd, Jffw = —b, and 

a = b + i(J?lw)o = -(Id + iH)J^w + i(je±w) . 

Define w = (Id + iH)w. It follows that w is holomorphic and odd, and 

3HPw = (Id + iH)J^w - i[H, JT]w = (Id + iH)J^w - i(J^±w) 

= —a 

as required. □ 
We will study the operator appearing in Lemma 14,41 

5 : C™(d + S(M)) -)■ C°°(M), Sh = (J^h^o- (4.1) 
The main point is the following result. 

Lemma 4.5. The operator S : Cg>(d+S(M)) -> C°°(M) is surjective. 

Given this, it is easy to see that any attenuation coefficient admits holo- 
morphic and antiholomorphic integrating factors. 

Proof of Proposition \4-l\ Given the attenuation a, we can take b to be any 
antiholomorphic even function in C°°(SM) (for instance 6 = 0). Choose 
h £ C™{d + S(M)) such that iSh = a-b - i(Jfiu b ) . Then 

b + i(J4?j_u b )o + i{Ji?±_h$)a = a 

and Lemma 14.41 shows that M'w = — a for some odd holomorphic w. The 
antiholomorphic case is analogous. □ 

It remains to prove Lemma 14.51 We will establish the surjectivity of S 
by proving that its adjoint is injective and has closed range. In fact, the 
adjoint involves the unattenuated ray transform 1° on 1-forms. Recall the 
space L 2 (dj r S(M)) where fJ<(x,£) = —(£,v(x)) and 

(h,ti) LUa+sm) = [ hh'fxd(dS(M)), 
M Jd+S(M) 

and d(dS(M)) is the natural volume form on dS(M) [6j Appendix A. 4]. 
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Lemma 4.6. If h £ C™(d + S(M)), f £ C 6 °°(M int ) then 

{Sh,f) L 2( M ) = (h, -— I°^1/)l2(s + 5(m))- 
Proof. We claim that for any v € C°°(SM) one has 



(Jf L f)vd(SM) = - f{J>t? ± v)d(SM). (4.2) 

SM JSM 
If this holds then we obtain by Santalo's formula [6, Appendix A. 4] that 

[(Sh)fdM = ^-[ (je^)fd(SM) = -^[ h^(Jif ± f)d(SM) 
Jm 27T Jsm 27r JSM 

= ~^f . h^ t (x,0)^±f&t(x,t))vdtd(dS(M)) 

1 



d+S(M) JO 

h(x,Ol°^f(x,Ovd(dS(M)) 

d+S(M) 



which is the required result. 

To prove (j4.2[) it is enough to show that 

I r ^±vd(SM) = 
Jsm 

for any v £ C°°{SM) with u = near <9M. Using the isometry F : SM — > 
SM, h-> (a;, — and the invariance of V, a change of variables and 

Santalo's formula imply that 

Jf ± vd(3M)= f (^yVjvdiSM) 
sm Jsm 

[ £ j F*(Vjv)d(SM) = [ t j Vj(F*v)d(SM) 
Jsm " Jsm 

Jf(F*v) d(SM) 



SM 

t(*,0 



d+S{M) JO 



Jf(F*v)(<p t (x, d* d(dS(M)). 



Since Jf?(F*v)((p t (a:, £)) = ^[(F*u)(<^(a;,£))] and F*v vanishes near dM, 
the last integral vanishes. □ 

Note that J>t±f(x,£) = = *df{£), where we write a(£) = oyf' 

for a 1-form a and tangent vector £, and * is the Hodge star operator. Thus 
the formal adjoint of S is given by 

S* : C c °°(M int ) C™(d+S(M)), S*f = --*-/ (*#(<£))■ 

Z7T 

The injectivity of S* follows immediately: if S*f = then the ray transform 
of the solenoidal 1-form *df vanishes, which implies that / = by pQ. To 
prove surjectivity of S it would then be enough to show that S* has closed 
range in proper spaces. The actual proof of the surjectivity will be slightly 
different, and we will proceed as in [231 Theorem 4.3]. 
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Proof of Lemma \4-5\ We may assume that (M, g) is embedded in a com- 
pact surface (N, g) without boundary, and that there is a finite open cover 
{U u . . . , U k } of N such that M C Ui, M HUj = for j > 2, and each 
(Uj,g) is simple. Let ipj G C^°(Uj) be a partition of unity so that tpj > 0, 
X^i=i = 1 in iV, and 991 = 1 near M. Let also Ij be the geodesic ray 
transform (with zero attenuation) on 1-forms in (Uj,g). Define the operator 
acting on smooth 1-forms on N, 

k 

P : C^iN^ 1 ) -> C^iV.A 1 ), Pa =Y t <p j i;i j (<p j a). 

i=i 

For the following details see [231 Theorem 4.3]. The principal symbol of P 
is given by the following expression, where C2 is a constant, 

If —A is the Laplace-Beltrami operator on N mapping smooth functions to 
smooth functions whose integral over N vanishes, define the operator 

A : C°°(N, A 1 ) ->■ C°°(N, A 1 ), Aa = -c 2 d(-A)-^ 2 5a. 

Then P + A is an elliptic pseudodifferential operator of order — 1 acting on 
1-forms in N. 

Consider the operator 

T : C°°(N,A 1 ) C°°(M), Ta = r M Sih(tp l a). (4.3) 

Here ru is the restriction to M, and Si is the operator (|4.ip in (Ui,g). We 
wish to show that T is surjective. Assuming this we can easily finish the 
proof of the lemma: If / G C°°(M) is given, there is a smooth 1-form a in 
N with tmSiIi (<pia) = f. Define 

h( x ,o = (JifaicO^foO, (x,o e a+,s(M), 

where uy, x is the map (I2.ip in (Ui,g). Since M is strictly contained in Ui 
we have h G C%f(d+S(M)), and it follows that Sh = tmSiIi^io) = / as 
required. 

It remains to prove surjectivity of (j4.3|) . In order to do this we express T as 
the adjoint of an operator involving P. The dual of C°° [M ) may be identified 
with the set S>' M {N) = {v £ 3>'(N) ; supp(v) C M}. Given v G @' M (N) 
choose Uj G C°°(N) with — ► u in @'(N) and supp('Uj) C {ipi = 1}. 
Then for a G C°°(A r , A 1 ), in the dual pairing in the indicated manifolds, by 
Lemma 14.61 we have 

(Ta,v) M = lim ((piSiIi((pia) > Vj) 17i 

= ~^ lim ( I i(Vl a )> I i(* d ( ( Pi v j)))lJf i (d + S(u 1 )) 

1 1 
= -—lim (a,P(*dvj)) N = - — (a,P(*dv)) N 

since P is continuous on @'(N). Now A(*cfo) = 0, so we consider 
Q : @'(N) -)■ ^'(^A 1 ), = __(p + A)(*d«). 
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It follows that T = r M Q\ since for a G C°°(N, A 1 ) and v G C c °°(M int ) 
(Toc,v)m = (a,Qv) N = (Q*a,v) N = (r M Q*a,v) M - 

Thus T in (|4.3p is a continuous linear map between Frechet spaces. Its 
adjoint is given by 

T* = Q\ 3 , M{N) :®' M {N)^@\N,k l ). 

The operator T* is injective. To see this, let v G @>' M {N) satisfy T*t> = 0. 
Then also (P + A)(*dv) = 0, showing that dv is smooth by elliptic regularity 
for P + A. Thus v is smooth and (^1-^1-^1 (^i * dv) = 0. Consequently 
IiifPl * cfo>) = in d+S(Ui), and the uniqueness result for I\ [I] implies 
that = dp near M for some smooth function p. Since v is smooth and 
supported in M, this is only possible if v = 0. 

Using that Q is an elliptic pseudodifferential operator, a standard argu- 
ment (see for instance [SJ proof of Theorem 6.3.1]) shows that the range of 
T* is weak* closed in ^'(N, A 1 ). It follows that T is a continuous linear 
map between Frechet spaces, and T* is injective with weak* closed range. 
Then T is surjective by [29\ Theorem 37.2], which concludes the proof. □ 

5. HOLOMORPHIC SOLUTIONS AND UNIQUENESS 

In this section we will prove the uniqueness result establishing injectivity 
of the attenuated ray transform. By using Proposition 14. H the boundary 
problem 

(Jf? + a)u = -f in SM, u\ 9S (m) = 
is equivalent with the problem 

J4?( e - W u) = -e~ w f in SM, e- w u\ ds{M] = 0. 

Here w is holomorphic, so also the right hand side e~ w f is holomorphic. 

The next main ingredient in the uniqueness proof is the following result, 
showing that in this situation also the solution e~ w u is necessarily holomor- 
phic. 

Proposition 5.1. Let (M, g) be a 2D simple manifold, and let / be a 
smooth function on SM holomorphic (resp. antiholomorphic) in the angular 
variable. Suppose that v is a smooth function on SM satisfying 

Jifv = -f in SM, v\ ds(M) = 0. 

Then v is holomorphic (resp. antiholomorphic) in the angular variable, and 
^o = 0. 

Proof. We only consider the holomorphic case. It is enough to show that 
(Id — iH)v = 0. The commutator identity (|2.3p shows that 

(Id - iH)v = (Id - iH)jev + i[H, J?]v 

= -(Id-iH)f + iJe ± v +i(J?±v) . 

We have (Id — iH)f = fo since / is holomorphic. Now 

(Id - iH)v = -h- (5.1) 
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where h = fo — i(J%±v)o is a function and a = *d(—ivo) is a solenoidal 
1-form (* is the Hodge star operator). 

What is important here is that there is no ^-dependence in h and a. In 
fact, since (Id — iH)v vanishes on dS(M), (|5.ip implies that the unatten- 
uated ray transform of h + cij£? vanishes identically. Uniqueness in the 
unattenuated case [T], using that a is solenoidal, implies that h = and 
a = 0. We have proved that 

J4?(Id-iH)v = 0, 

and (Id — iH)v = since v\qs(m) = 0- d 

It is now easy to give a proof of the injectivity result for the attenuated 
ray transform of functions. 

Proof of Theorem{Tj\ Let first / € C c °°(M int ), and assume that I a f = 0. 
Then the function 

u(x,0 = r ( "°/(7(t,x,0)e /ota(7(s ' x ' ?))ds dt 
Jo 

is smooth in SM and satisfies 

(JT + a)u = -f in SM, u\ ds{M) = 0. 

By Proposition 14.11 there is a holomorphic w € C°°(SM) and an antiholo- 
morphic w G C°°(SM) such that 

jr(e-™u) = -e~ w f, 

J?(e-™u) = -e~™f. 

Now Proposition 15.11 shows that e~ w u is holomorphic and e~ w u is an- 
tiholomorphic. Multiplying by e w and e w , it follows that -u itself is both 
holomorphic and antiholomorphic. This is only possible when u does not 
depend on £, that is, u = uq. Now the transport equation reads 

du (0 + au = -/ in SM, ilo\qs(m) = 0. 

Evaluating this at ±£ and substracting the resulting expressions gives that 
duo = 0, hence no = by the boundary condition. Consequently / = 0. 

It remains to prove the result when / € C°°(M) may have support ex- 
tending up to the boundary. In fact, this case can be reduced to the result for 
compactly supported functions by using the general principle that (I a )*I a 
is an elliptic pseudodifferential operator. 

Suppose / € C°°(M) and I a f = 0. We consider more generally the 
weighted ray transform with weight p E C°°(SM), 

Ipf(x, = [ T{X '° P(<Pt(x, 0)/(7(*, ^, 0) dt, (x, G d + S(M). 
Jo 

With the choice p = e~ u - , we obtain I p f = 0. Let (M,g) 3D (M,g) be 
another simple manifold which is so small that any M-geodesic starting at 
a point of d-S(M) never enters M again. We extend a to M as a smooth 
function and / by zero to M, and denote by I p the corresponding weighted 
ray transform in M. 
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It is easy to see that I p f(y, rj) = for all (y, 77) G d + S(M), since either the 
geodesic starting from (y, 77) never touches M or else I p f(y, rj) = I p f(x, £) for 
some (x,£) € 5 + 5(M). By [TUl Proposition 2], since /j is nonvanishing, I*I p 
is an elliptic pseudodifferential operator of order —1 in M . Now I p I p f = 0, 
and elliptic regularity shows that / is smooth. Thus / G C^°(M mt ) and 
I p f = 0, showing that I a f = 0. The result above implies that / = as 
required. □ 

We proceed to the attenuated ray transform of 1-forms. In this case the 
required generalization of Proposition 15.11 is as follows. 

Proposition 5.2. Let (M,g) be a 2D simple manifold and let p G C°°(SM) 
be holomorphic (resp. antiholomorphic) in the angular variable. Suppose 
F(x,£) = f(x) + aj(x)£ 3 where / is a smooth function and a is a smooth 
1-forrn on M. If v G C°°(5'M) satisfies 

JTv = -pF in SM, v\ ds(M) = 0, 

then v is holomorphic (resp. antiholomorphic) in the angular variable. 

Proof. We only prove the holomorphic case. As in Proposition 15. 11 we study 
the function (Id — iH)v and note that 

JT(Id - iH)v = (Id - iH) J^v + i[H, J>t?]v 

= -(Id - iH)(pF) + iJ^ ± v + i(Jfiv)o. (5.2) 

We would like to show that (Id — iH){pF) is a first order polynomial in £. 
The first step is to prove that 

(Id - iH)((p - Po )F) = ~h (5.3) 

for some h = h(x). This can be reduced to elementary facts about Fourier 
coefficients. Fix a point x in M, and choose an orientation preserving isom- 
etry $ : S x — > S 1 . If Hq is the Hilbert transform on S 1 , it is enough to prove 
the equivalent statement 

(Id - iHo)^- 1 )*^ - po)F) = h(x). (5.4) 

But (Id — iHq)^^^ 1 )* (p — po)) = since p is holomorphic, meaning that 
one has the Fourier series 

00 

^- i r( P -po) = Y J ^ ke . 

k=l 

Also, since F is a first order polynomial in £, one has the Fourier series 
((J)- 1 )*^ = 6 + bie 10 + b_ie~ td . Multiplying these Fourier series gives (|5.4p . 
By (pT2|) and (fOj) . we obtain 

J^(ld - iH)v = -p (Id - iH)F - F 

where F = h — i(J^t>)o — i * dfo(0 is a first order polynomial in £. Next 
we employ a Hodge decomposition a = dp + *dq where p,q G C°°(M). One 
has the antiholomorphic projections 

(Id - iH)(dp(£)) = (Id - iH)JTp = Jf(Id - iH)p - i[H, Jf]p 
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and similarly 

(Id - iH)(*dq(£)) = (Id - iH)J^ ± q = (Id - iH)[H, J?]q = (Id - iH)HJ?q 
= HJ^q + i(J^q - (M 3 q) ) = i(Id - iH)J?q 
= i(je-iJ^ ± )q, 

the last line using the above computation for p. 

Putting these results together, we have proved that 

JT(Id - iH)v = —F in SM, (Id - iH)v\ dS{M) = 0, 

where F is the first order polynomial in £ given by 

F = pof + h- i(J#±v) + po(d -i*d){p + iq)(£) - i * dv (£). 

This shows that the geodesic ray transform of F vanishes. Therefore F(£) = 
dp{^) for some p € C°°(M) with p\qm = 0, and the equation for (Id — iH)v 
implies that 

(Id — iH)v = —p. 

This proves that v is holomorphic. □ 

The injectivity result for 1-forms is now proved in a similar way as the 
corresponding result for functions. 

Proof of Theorem \1.2l Let / be a smooth function and a a smooth 1-form 
on M, and let F = f + ayf. Assume first that / and a are compactly 
supported in M mt . Then the function 

u(x,i) = ! T{X * F{i Pt (x,0)e f ° a{ ' ris ' x ' 0)ds dt 
Jo 

is smooth in SM and satisfies the equation 

(JP + a)u = -F in SM, u\ ds{M) = 0. 

As before we use Proposition 14.11 to find a holomorphic function w and an 
antiholomorphic function w such that 

J?( e - W u) = -e~ w F, 

^(e-^u) = -e-*F. 

Proposition 15.21 implies that e~ w u is holomorphic and e~ w u is antiholomor- 
phic, and consequently u is both holomorphic and antiholomorphic. Thus 
u = uq, and the transport equation can be written as 

(duo + + au = -f in SM, ii |as(M) = 0. 

Consequently a = —duo and / = —auo which is the required result. 

Again the general case where / € C°°(M) and a is a smooth 1-form in 
M may be reduced to the previous case by elliptic regularity. Assume that 
I a (f + = 0. If p = e~ u - , this implies that I p (f + ctj&) = where 

I p F(x,0 = [ ( " { p(<pt(x,0)F(Mx,0)dt. 
Jo 
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Consider the solenoidal decomposition a = a s + dp, where 5a s = and 
— Ap = 5a with p\qm = 0. An integration by parts shows that we have 
l p (f-ap + a^) = 0. 

Let (M,g) DD (M,g) be a simple manifold as in the proof of Theorem 
II. 1| and extend a smoothly to M and f,p,a s by zero to M. It follows 
that I p (f — ap + a s ^ 3 ) = 0. By [TTJ Proposition 1] (where we make the 
choices w = p and a = p, so that the modified elliptic condition in [TTJ 
Remark 1] is satisfied), I p Ip is a pseudodifferential operator of order —1 in 
M mt which is elliptic in the sense that whenever f',a' are in L 2 (M mt ) and 
I*I p (f '+a'j(, 3 ) = and 5a' = 0, then /' and a' are smooth. Thus, f — ap and 
a s are smooth and compactly supported in M mt and I p (f — ap + a s ^ 3 ) = 0) 
showing that f — ap = ap and a s = dp for some smooth p with p\ d ^j = 0. 
Since a s is zero outside M it follows that p vanishes outside M, and one 
obtains / = a(p+p) and a = d(p+p) in M with p+p vanishing on dM. □ 

Finally we give a proof of the stability result in the introduction. 

Proof of Theorem \1.'J[ We make the same preparations as in the end of proof 
of Theorem 11.21 Thus, we consider the decomposition a = a s + dp with 
p = G{5a) in M, choose a slightly larger manifold (M, g), and extend a 
to M and f,p,a s by zero to M. Let iV = I*I p with p = e~ w - . Since the 
modified elliptic condition of [IT] Remark 1] is satisfied, [TTJ Proposition 2] 
implies the estimate 

11/ - ap\\ L 2 {M) + ||a s || L 2 (M) < C{\\N(f -ap + a^lli^M) 

+ 11/ - a P\\ H -HM) + IKHtf-i(M))- ( 5 - 5 ) 

Let L 2 S (M) = {cp+Pjtf ; <j> G L 2 (M), (3 is a 1-form in L 2 (M) and 5(3 = 0} 
be the space of solenoidal pairs, assumed to be extended by zero to M. As 
in the end of proof of Theorem 11.21 one can use the ellipticity of and 
injectivity of I p to see that N : L 2 S (M) — > H l {M) is a bounded injective 
operator. We can then use |11| Lemma 1] with the choices X = L 2 (M), 
Y = H l (M), and Z = F _1 (M) (the latter two being the natural Sobolev 
spaces for solenoidal pairs), and conclude from (|5.5|1 that 

11/ - ttp||x2(M) + IKIlL 2 (Af) < C\\N '(/ - ap + a^)[| H i(M)- 

The stability result follows by noting that N(f + Ojf 7 ) = N(f — ap + a s ^ 3 ) 
and by taking M X =M. □ 

6. Reconstruction procedure 

Let (M, g) be a simple 2D manifold and let a be a smooth complex func- 
tion on M. In this section we give a procedure for determining a smooth 
function / in M from the knowledge of I a f. 

There are two nontrivial parts in the procedure: computing the inverse 
of the unattenuated ray transform 1° in (M, g), and the construction of 
(anti)holomorphic integrating factors for the equation (Jif + a)u = 0. If 
(M, g) has constant curvature these operations can be done explicitly (since 
the W operator vanishes, see [23] and Corollary I4.3[) . Also, if (M,g) is a 
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small perturbation of a constant curvature manifold then the W operator has 
small norm on L 2 (M) and these two operations can be expressed in terms 
of convergent Neumann series, see |14j and Corollary 14.31 again. However, 
for general simple (M, g) it is not clear how to carry out these operations in 
an explicit way. 

For simplicity we will assume below that / € C^°(M int ), since in this case 
all the functions will be smooth up to dM and we do not need to worry about 
regularity issues. Theorem II .41 immediately follows from the next result. 

Proposition 6.1. Under the stated assumptions, a function / € C^°(M mt ) 
can be determined from the knowledge of I a f using the following procedure: 

1. Define a function d on dS{M) by 

d(x t ) = i l a f(*,0, M)€d + S(M), 
I 0, otherwise. 

2. Find a holomorphic function w and an antiholomorphic function w, 
both smooth odd functions on SM, such that Jtfw = Jtfw = —a. 

3. Let = (Id - iH){e~ w d) and /3 = (Id + iH)(e~™d) on dS{M). 

4. Let v = f3 o i/) + u^W* and ?) = /3 o ^ + u {l0) ~ 1{A1 ^ in SM, 
where A*_f3 = (3 — (3 o ifi on d±S(M) and (I ) -1 is the inverse of the 
geodesic ray transform in (M,g), in the sense that 

(I°)- 1 I°(<l> + a j ?) = <l> + a j ? 

for a smooth function <p and a solenoidal 1-form a. 

5. Define rh = |(Id — iH)(e w v) + |(Id + iH)(e w v) and u = rh — mo- 

6. Define q = (d — u)q o ip + (u^" +aM )-)o, and let u = q + u. 

7. Let / = — (Jfu + au)o. 

Proof. Let u be the solution of J$?u + au = —f in SM with u|,9_s(m) = 0, 
so that u\d + s(M) = I a f ■ If w an d ?5 are as described, then one has the two 
equations 

je(e- w u) = -e~ w f, 
je(e-™u) = -e~™f. 

The right hand side in the first equation is holomorphic in the angular 
variable. We will show that v = (Id — iH)(e~ w u), the antiholomorphic part 
of the solution e~ w u, is determined by I a f. In fact, the computation in 
Proposition 15. II shows that 

Jtfv = —cj) — Olj^ 3 

where <f> = f — i(J4?±(e~ w u))o and a = *d(—i(e~ w u)o). Here we used that 
( e ~ w )o = 1 since w is odd. We have that v — v o ip\g + s(M) = ^ (<^ + a j&)- 
Here a is solenoidal, so 1° is invertible and 

(f> + a^ j = (I°)-\v-voil>\ d+s{M) ). 

This proves that v is the function given in Step 4 above. A similar argument 
shows that v = (Id + iH)(e~ w u) is the other function in Step 4. 
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We have obtained two decompositions 

e~ w u = h + —v, 
e~ w u = h + —v 

where h is holomorphic, h is antiholomorphic, and v and v can be determined 
from the attenuated ray transform of /. This results in two decompositions 
for the solution u, 

u = e ™h + ^e w v, 
u = e ih h + ^e™v, 

where again e w h is holomorphic, e w h is antiholomorphic, and e w v and e w v 
are known. This determines u up to a term which is constant in £, which 
can be seen by writing 

u = l(Id + iH)u + ^(Id - iH)u = ^(e w h) + ^(e*h) + m 

with m given in Step 5. 

Write u = q + u where uq = 0. Then q = uq and u = m — mo- To find the 
term q we note that g|aM = u — u\qm = (d — u)q\qm, where necessarily u — u 
is independent of £. Taking the odd part in the equation it + au = — / 
implies that 

jTg = -(JTii + an)_. 

Therefore q is given by the quantity in Step 6. We have determined the 
solution u in SM from the knowledge of I a f. Now / = —Jtfu — au, and 
taking averages proves the formula in Step 7. □ 
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